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Abstract Wc investigate the background dynamics 
when dark energy is coupled to dark matter in the uni- 
verse described by loop quantum cosmology. We con- 
sider dark energy of the form modified Chaplygin gas. 
The dynamical system of equations is solved numeri- 
cally and a stable scaling solution is obtained. It hence- 
forth resolves the famous cosmic coincidence problem in 
modern cosmology. The statefinder parameters are also 
calculated to classify this dark energy model. 



1 Introduction 

Recent observations of type la Supernov ae indicate that 
Universe is expanding with acceleration ( Perlmutter et al^ 
19991 Riess et al 1998 ) and lead to the search for a new 
type of matter which violates the strong energy con- 
dition, i.e., p + 3p < 0. In Einstein's general relativ- 
ity, an energy component with large negative pressure 
has to be introduced in the total energy density of 
the Universe in order to explain this cosmic acceler- 
ation . This energy component is known as dark en - 
ergy (jSahni &: StarobinskviuOOOt iPadmanabhanI 120031 ) . 
The re are many candid ates supporting this behav- 
ior (ICopeland et a scalar field or quintessence 

( Peebles fc Ratralbood ) being one of the most favored 
candidates as it has a decaying potential term which 
dominates over the kinetic term thus generating enough 
pressure to drive acceleration. 
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Presently we live in an epoch where the densities 
of the dark energy and the dark matter are com- 
parable. It becomes difficult to solve this coinci- 
dence problem without a suitable interaction. Gen- 
erally interacting dark energy models are studied to 
expla i n the cosmic coincidence problem ( Jamil et~all 
2OIOI: Ijamil fc SaridakisI boiot Ijamil fc Faroool boiot 



Jamil et all I2OIOI : I Jamil k RahmanI I2OO9I) . Also the 
transition from matter domination to dark energy dom- 
ination can be explained through an appropriate energy 
exchange rate. Therefore, to obtain a suitable evolution 
of the Universe an interaction is assumed and the decay 
rate should be proportional to the present value of the 
Hubble parameter for good fit to the expansion history 
of the Univ erse as determined by the Su pernovae and 
CMB data Jjamil et allboid Ijamil fc Saridakis,.201ol) . 
A variety of interacting dark energy mod els have been 
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proposed and studied f or this purpose (ISetar 
Hu fc Lind hood IWti fc YtJ l2007i Ijami 



2007 



2010; 



In recent years, the model of interacting dark energy 
has been explored in the framework of loop quantum 
cosmology (LQC) as well: It is shown in (|Wu fc Yu 
20081) that for the quintessence model, the cosmological 
evolution in LQC is the same as that in classical Ein- 
stein cosmology, whereas for the phantom dark energy 
the loop quantum effect significantly redu ce the pa- 
rame ter spacetime required by stability. In ( Chen et all 
20081 ) , the authors used a more general interaction term 
to study the interacting dark energy. They showed that 
in LQC, the parameter space for the existence of the ac- 
celerated scaling attractor is found to be sma her then 
that in Einstein cosmology. In another study (|Fu et al 
20081 ). the authors studied the model with an interact- 
ing phantom scalar field with an exponential potential 
and deduced that the future singularity appearing in 
the standard FRW cosmology can be avoided by loop 
quantum effects. 
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111 this paper, we extend the model of interacting 
modified Chaplygin gas (MCG) from the framework of 
Einstein gravity to LQC. We construct a dynamical 
system of equations and solve them numerically. We 
obtain a stable scaling solution (which is also an 'at- 
tractor') of modified FRW equations. We discuss our 
results in the final section. 



2 The model 



The modified Fr i edmann equation for LQC is gi ven by 
( Wu fc Yulliooi: IChen et aJliooj IFu et alllioosh 



3 V pi 



(1) 



Here pi = \/3n^^^G^h is the critical loop quantum den- 
sity and 7 is the dimensionless Barbero-Immirzi param- 
eter. We assume the interaction between dark energy 
and pressureless dark matter. Hence the energy bal- 
ance equations for the interacting dark energy and dark 
matter can be expressed as 



Pmcg + 'iH{l 



:g)Pmcj 

3Hp, 



(2) 
(3) 



where Q = 'ibHp is the interaction term, b is the 
coupling parameter (or transfer strength) and p — 
Pmcg+Pm is the total cosmic energy density which satis- 
fies /0 -H3g(/0-Hio) = |Guo fc ZhandboOsHCampo et all 
2008[) . Note that addition of the above two equations 
leads to the energy conservation. Due to unknown na- 
ture of both dark energy and dark matter, the interac- 
tion term can not be derived from the first principles. It 
is worthy to note that if Q < than it will yield the en- 
ergy density of dark energy to be negative at sufficiently 
early times, consequently the second law of the rmody- 
namics can be violated (jAlcaniz fc Limal 120051 ) hence 
Q must be positive and small. From the observational 
data of 182 Gold type la supernova samples, CMB data 
from the three year WMAP survey and the baryonic 
acoustic oscillations from the Sloan Digital Sky Survey, 
it is estimated that the coupling parameter between 
dark matter and dark energy must be a small positive 
value (of the order unity), which satisfies the require- 
ment for solving the cosmic coincid ence problem an d 
the second law of thermodynamics ([Feng et a 
Because of the underlying interaction, the beginning of 
the accelerated expansion is shifted to higher rcdshifts. 

Consequently we obtain the modified Raychaudhuri 
equation 



H = ~lip + p)(l-2^ 
2 V pi 



(4) 



where p is the total pressure {p = Pmcg)- We shall use 
modified Chaplygin gas as the dark energy. The MCG 
equation of state is given by 
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(5) 



where A, B and a are constants. The MCG best 
fits with the 3-year WMAP and the SDSS data with 
the choice of p arameters A = —0.085 and a = 1.724 
( Lu et aJ 2008 ) which are improved con straints than 
the p revious ones —0.35 < A < 0.025 (jjun fc Zhou 
20051 ). Recently it is shown that the dynamical attrac- 
tor for the MCG exists at Wmeg = —1, hence MCG 
crosses this value from either side Wmcg > — 1 or a;,ncg < 
— 1, independen t to the choice of model parameters 
( Jing et a I2OO8I). A generalization of MCG is suggested 
in |Debnathll2007t) by considering B = B(a) = Bou'', 
where n and Bo are constants. The MCG is the gcner- 
alization of generalized Chaplygin gas Pmcs = —B / p^cg 
( Barreiro fc SCTi |2004 ICarturan fc Finellil boOSf ) with 
the addition of a barotropic term. This special form 
also appears to be consistent with the WMAP 5— year 
data and henceforth the support the unified model with 
dark en ergy and matter based on generalized Chaply- 
gin gas jBarriero et allboOSt iMakler et a3l2003l: ISetare 



200a 2007 



a2007l ). In the cosmological context, the 



Chaplygin gas was first suggested as an alternative 
to quintessence and demonstrated an in creasing A be- 
havio r for the evolution of the universe (jKamenshchik 
20011 ). Recent supernovae data also favors the two- 



fluid cosmological model with Chaplygin gas and mat- 
ter (|PanotODOulosll2008l ). 

To analyze the dynamical system, we convert the 
physical parameters into dimensionless form as 



X = In a, u 



3ff2 



(6) 



where ao = 1 is assumed, where the subscript refers 
to the present time. Making use of ([T]) to dl]), we can 
write 



du 
dx 



dv 



-36(u + u) — 3u(l -t- Wmcg) 



+3u[u(l-t-Wmeg) +!'](- l + (7) 



— = 3b(u + v) — 3v + 3v[u(l + U!a^c<y) + 

dx ° 



X -1 



U + V 



(8) 



where the state parameter of modified Chaplygin gas is 

Pmcg , B{U + i;)2("+l) 



A 



Pn 



p^+iu"+i(w-hu- l)"+i' 



(9) 
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For the mathematical simphcity, we work out a = 1 
only. The eritieal points of the above system are ob- 
tained by putting = = 4^ which yield 



Ulc 



Vie = 



1-b ^{-l + A{-l + b)){-l + brpi 

bpi 

-B+{-l + A{-l + b)){-l + b)pl 

{-l + A{-l + b)){-l + b)pi 



+ ^B{-l + A{-l + b)){-l + b) 
1 \/B 



(10) 



U2c = 



V2c 



V{-l+A{-l+b)){-l+brpi 

bpi 



B-{-l + A{-l + b)){-l + b)pl 



i-l + A{-l + b)){-l + b)pi 



+ ^B{-l + A{-l + b)){-l + b) 



(11) 



The two critical points correspond to the era dominated 
by dark matter and MCG type dark energy and exist 
for A > ji^. 

For the two critical points, the state parameter (9) 
of the interacting dark energy takes the form 



1,2 



which holds only when Uic + Vic ^ 1. 

We further check the stability of the dynamical sys- 
tem (Eqs. (7) and (8)) about the critical point. To 
do this, we linearize the governing equations about the 
critical point i.e. u — Uc + 6u and v — Vc + Sv, we obtain 



it 



-3{-l + b + 2u + v) 

.4|-3 + 6uf-H 

3B{u + v)^ 
u^(—l + u + vYp\ 



[u + v) 



X ^2u^ - (-1 + v)v'^ - uv(2 + v) 

+ w)(-3 + 2v) + u^{~5 + 4w)| 
2Au 



3b~3v 1 



A{-2v + (u + v)"^)' 



{u + vy 



3Bv{u + vy' 



u^i-l + u + vypjl 
-(-2 + + v)v - u(-4 + f2)| 



5u 



35 + 



+ U + vy{u + vYpI 

-t-i;(8 + 3(-3 + v)v) + u{2 + w(-12 + 7v))) 
-u{-l + U + vf{u^{-l + A{-2 + u) + u) 
-t-2u(-l + (2 + A)u)v + (-1 + (5 + A)u)v^ 



+2v- 



Sv. 



The subscript c refers to quantities evaluated at the 
critical point of the dynamical system. We also calcu- 
late the deceleration parameter q = —1 — (H/H^), in 
this model as 



1 I " ( 1 I Pnic, 

2V p 



Pmcg \ / I - 2p/pi \ 
p p/pi ) 



(12) 



which can be written in terms of dimensionless density 
parameter fimcg = Pmcg/p: 



9LQC ^ + 2 ("^ ~^ ^^mcgf^ii 



Clearly in the limit of p\ — > 
for the Einstein's gravity as 



l-2p/pA 
l-p/Pi I' 

oo, we retrieve the result 



<]EG — ^1+2(1+ ^mcg^^mcg 



Assuming p/ pi ~ e ^ 0{1) and using ([B]), we obtain 

■1 - 2e 



..In 



meg 



U + V 



1 - e 



(13) 



Since the only physically acceptable solution corre- 
sponds to the first stable critical point, such that 
{u,v) (uicVic). Hence ([T^ gives 



-1 + -X, X 



/^^^w^ui^y (14) 



-35-3u 1 



{u + v) 



As special cases, observe that for e = 1/2, we have 
q ~ —1 while e = 1 yields q — > —00 if oJmcg > — 1/wic- 
Moreover the Hubble parameter varies as 



6B{u + v) 



'l—^Iu^ + 2u\-l + v) H=—, 



u{—l + U + vYp\ 
+ {-2 + v)v + u{\ + (-1 + v)v)^ 



(15) 



5v, 



where we have ignored the integration constant. Inte- 
gration of (fT5|) yields 



a{t) = aot^x ^ 



(16) 
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which gives a power law form of the expansion. 

We also calculate the statefinder parameters. Sahni 
et al (jSahni et all boosl ) introduced a pair of cosmo- 
logical diagnostic pair {r, s} which they termed as 
Statefinder. The two parameters are dimensionless and 
are geometrical since they are derived from the cosmic 
scale factor alone, though one can rewrite them in terms 
of the parameters of dark energy and matter. Addition- 
ally, the pair gives information about dark energy in a 
model independent way i.e. it categorizes dark energy 
in the context of background geometry only which is 
not dependent on the theory of gravity. Hence geomet- 
rical variables are universal. Also this pair generalizes 
the well-known geometrical parameters like the Hubble 
parameter and the deceleration parameter. This pair 
is algebraically related to the equation of state of dark 
energy and its first time derivative. 



on forever with an ever increasing rate. The variations 
of 9LQC, ''LQC and slqc are shown in figure 5-7 re- 
spectively against uj^cg- It is observed that the more 
negative the state-parameter of MCG, the more nega- 
tive values will be taken by the deceleration parameter. 
Finally our results also reduce to those for Einstein's 
gravity under suitable limits of parameters. 
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am' ' 3{q-l/2)' 
In the present model, (fT7)) gives 

^LQC = (^1-^)(1-3X), 
SLQC = 2X. 



(17) 



(18) 
(19) 



It is interesting to note that the pair {^lqc, slqc} 
yields the ACDM (cosmological constant-cold dark 
matter model) {rEcSEc} = {IjO} when X = (or 
e=l/2). 



3 Discussion 



In this work, we considered modified Friedmann model 
in loop quantum cosmology. We assumed dark energy 
of the form modified Chaplygin gas. The interaction 
between dark matter and MCG has been investigated 
in LQC. The dynamical system of equations is solved 
numerically and a stable scaling solution is obtained. 
It henceforth resolves the famous cosmic coincidence 
problem in modern cosmology. The deceleration pa- 
rameter and statefinder parameters are also calculated 
to classify this dark energy model. The dimcnsionally 
density parameters v and u are drawn in figures 1 and 
2. We see that v decreases and u increases during evolu- 
tion of the universe. From figure 4, we also see that the 
ratio of the above parameters decreases during time. 
The phase space diagram (figure 3) shows the attractor 
solution hence the present state and the future evo- 
lution of the universe is independent to the choice of 
initial conditions. Moreover the expansion of the uni- 
verse is governed by a power-law form, rather than ex- 
ponential or oscillatory. Hence the expansion will go 
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Fig. 1 The dimensionless density parameters are plotted 
against e-folding time. The initial condition is v{0) — 0.9, 
u{0) — 0.2. Other parameters are fixed at 6 = 0.2, A = 0.3, 
B = 0.5 and a = 1. 




Fig. 2 The dimensionless density parameters are plotted 
against e-folding time. The initial condition is v{0) = 0.6, 
m(0) = 0.6. Other parameters are fixed at & = 0.03, A = 0.3, 
B = 0.5 and a = 0.5. 
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6 



0.8- 



0.7- 



0.6- 



0.5- 



0.4- 



0.3- 



\ I / / 

-^^^^"Os>^ \ \ \ I y ^ 
^^"^H^N \ \ / 

^ ^ ^ "^?«^^ ■^Nw\ I / 

^5*^^ --^ ^ ^^S^X 
^ "^^1^^ ^ ^ ^^^>^ 
"-^^i^^ ^ ^ ^ ^ 

X ^ ^ ^ ^ ~ 

^ ^ ^ ^ ^ " 

N X X ^ 



0.6 



— I — 
0.7 



0.8 



0.9 



1.0 




Fig. 3 The phase space diagram of parameters depict- 
ing an attractor solution. The initial conditions chosen are 
i;(0) = 0.5, 'u(O) = 0.6 (green); d(0) = 0.6, u(0) = 0.6 
(blue); «(0) = 0.7, •u(O) = 0.6 (red); i;(0) = 0.8, m(0) = 0.6 
(brown). Other parameters are fixed at 6 = 0.3, A — 0.3, 
B = 0.5 and a = I. 



Fig. 5 The deceleration parameter is plotted against the 
state parameter. Other parameters are fixed at & = 0.3, 
e = 0.4, A = 0.4 and B = 0.5. 





Fig. 4 The ratio of density parameters is shown against 
e-folding time.. The initial condition chosen is v{0) = 0.6, 
u{0) = 0.6. Other parameters are fixed at 6 = 0.03, A = 0.3, 
B = 0.5 and a = 0.5. 



Fig. 6 The statefinder parameter r is plotted against the 
state parameter. Other parameters are fixed at b = 0.3, 
e = 0.4, A = 0.4 and B = 0.5. 
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Fig. 7 The statefinder parameter s is plotted against the 
state parameter. Other parameters are fixed at & = 0.3, 
e = 0.4, A = 0.4 and B = 0.5. 



